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ABSTRACT 


A class of nonlinear systems, which can be represented 
by a feedback system with one or two linear plants, and one 
single-valued time-invariant nonlinearity, is studied. The Volterra 
series describing such a system, where the nonlinearity can be 
represented as a power series, is shown to be unique. Initial 
conditions on the linear plants are incorporated into the Volterra 
series in a simple manner. Using this initial condition method, 
the convergence of the series is assured for inputs containing a 
step type function. The convergence of the series is investigated, 
and it is shown that for first and second order systems the series 
converges, and bounded-input bounded-output stability conditions are 
obtained for a nonlinearity whose slope lies in the Hurwitz sector. 
For higher order systems without zeros in their linear plants, it is 
shown that a contraction can be obtained by use of another norn, 
which also gives a bound on the output, for a nonlinearity whose 


slope lies in the Hurwitz sector. 
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CHAPTER I 


INTRODUCTION 


Nonlinear systems play an important part in the practical 
physical world. Every transducer or amplifier which is used to make 
up a system has some inherent nonlinearity especially when the signal 
it is called upon to handle reaches a sufficiently large magnitude. 
It would seem that systems could not be analyzed without considering 
the nonlinearities but fortunately the usual practice is to run the 
elements of the system such that they approximate linear devices very 
closely, Then it is possible to analyze and synthesize systems using 
linear devices recognizing the contraints on the signals which the devices 
will be called upon to process. However some systems, such as a phase 
locked loop Plies er a demodulator [29] base their operation upon a 
nonlinear device and in fact although designers of systems try to avoid 
using nonlinear devices due to the difficulty in estimating their effect, 
it may be that the use of such a device would enhance the properties of 
the system under design. There are also the cases of devices which are 
linear at small signals for which it is not practical to prevent the 
device from receiving a large enough signal to drive it into its non- 
linear region. A simple case of this type would be an audio amplifier 
te We Wate 

The ability to analyse and even synthesize a nonlinear system 
would be a valuable asset to the systems design engineer. This ideal is 


far from being reached due to the complexities involved in even the 


t+ Numbers placed within square brackets denote references. 
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simplest of nonlinear systems. 

A nonlinear system is defined as any system for which super- 
position does not hold. As such it can be a system with a time varying 
nonlinear element, a time invariant single valued nonlinear element, 
an element with hysteresis, or any combination of these types of non- 
linearity. Ideally, a theory for the analysis of nonlinear systems 
would be able to handle any nonlinearity, and the linear case would be 
just a special limiting case in this theory. This ideal too is far 
from fruition. 

There are two areas of interest to the systems engineer about 
the system under c@dstd deaendive The first one is whether the system is 
stable, either asymptotically stable or bounded-input bounded output stable. 
If the system is stable, then a measure of the performance of the system 
may be desired. For example, the rise time in response to a step input, 
or the steady state error of the system, may be important. Usually these 
two areas are treated by different methods, and there have been quite a few 
methods developed for different types of systems and different types of non- 
linearities. Gibson [12] in his book gives a good introduction to a 
variety of methods, most of which are limited to systems which can be 
_ described by first or second order nonlinear differential equations. 

The siete of asymptotic stability has been treated by 
Liapunov's second method (see [12]) and by Popov [see [1]). These 
methods give sufficient conditions for stability and in some special 
cases give necessary conditions too. For these methods the nonlinear 
function is single valued and contained in a sector [k, sk], and gives 


a zero output for a zero input (see [1]). 
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For nonlinear systems asymptotic stability does not guarantee 
stability with an input (see [12]) so the concept of bounded-input 
bounded-output stability is investigated. Although other methods such 
as describing function methods [12] have been used to investigate 
stability of systems with inputs, it appears that a functional analysis 
approach to this problem gives very general results quite easily. Also, 
using a functional series solution to the system equation, called the 
Volterra series, the system can be represented in terms of its parameters 
and the input function, for a large class of systems. 

V. Volterra [26] introduced the concept of a functional early in 
this century, and also introduced the type of integral equation with which 
most engineering systems are described. There have been many mathematicians 
working on this topic since then, but it was not until N. Wiener in 
1942 [27] applied the functional series to the study of a nonlinear elec- 
trical circuit problem that engineers saw the possibilities of the theory 
in the analysis of nonlinear systems. Since then Barrett [2], Brilliant 
[3], George [11], Zames [29], Parente [18], Christensen [4], Sandberg [22] 
and Holtzman [14] to mention just a few have investigated the stability of 
systems (to inputs) and considered the requirements for convergence of 
the Volterra series. In most cases it is found that the stability limits 
are very conservative compared to other methods although Zames [30], [31] 
has proposed a circle theorem which with some restrictions on the non- 
linearity gives a condition similar to that for the Popov criterion. 
However, a system which is bounded-input bounded-output stable may 
exhibit limit cycle oscillation when the input returns to zero. 

All these methods give only sufficient conditions for stability 


whereas in the linear systems necessary and sufficient conditions follow 
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from the application of the Nyquist or other similar criterion. This 
has led to the speculation that if the nonlinear function f(y) was 
replaced by a linear gain ky, and the range of k for which this linear 
system is stable [k, sk] (the Hurwitz sector) is determined, then the 
system may be stable for all nonlinear functions f(y) lying in the 
sector [k, ko] in the (y,f(y)) plane. This conjecture was originally 
put forward by M. Aizerman [1]. It has been shown to be true for a 
range of systems, but not true for many others. In fact Willems [28] 
has found that a counter example to the Aizerman conjecture can be 
found for almost any sufficiently smooth nonlinearity. However the 
Hurwitz sector gives an upper bound on the region over which a non- 
linearity may be expected to give a stable system. 

A more limiting criterion which will be used here is that the 
slope of the nonlinearity (where it exists) should lie within the range 
of [k, so], as defined by the Hurwitz sector for the system. 

| The second area of interest to the engineer, namely the per- 
formance of the system when subjected to different inputs, can be 
approximated by using the Volterra series. Also the response of the 
system to initial conditions can be approximated using the same method. 

In this thesis a class of systems will be examined using 
operator notation, where an operator will be signified by an underlined 
capital letter. The class of systems will be the class of feedback 
systems that can be represented by up to two linear plants with a single 
valued time invariant nonlinearity which passes through the origin, and 
which can be represented arbitrarily accurately by a finite power series. 
It will be shown with such a system that the Volterra series is a unique 


series expansion and that there is convergence of the series for non- 
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linearities whose slopes are within the Hurwitz sector. It will also 
be shown that for first and second order systems that bounded-input 
bounded-out put peices exists for such nonlinearities. A method of 
incorporating initial conditions into the series is also discussed, and 
using this method a means of getting better convergence of the series 
with step inputs is introduced. 

Chapter II introduces the notation and mathematics to be used. 
Operator algebra is discussed and the general system to be analysed is 
introduced. The space in which the inputs and outputs belong is then 
described, along with a weighted norm which is used in the convergence 
proofs. Some norms of operators are examined and the basic contraction 
mapping theorem is stated. 

Chapter III demonstrates the solution of the system equation 
using the Volterra series by two different methods,and then a proof is 
offered to show that these two methods give identical series in the 
limit. Then a simple way of evaluating the series using transform 
theory is developed for the general system under study followed by an 
example. 

In Chapter IV the convergence of the series is proved for 
functions belonging to the L., space. The first and second order systems 
are shown to give bounded-input bounded-output stability for nonlinear- 
ities as mentioned above and for higher order systems the output is shown 
to be exponentially bounded. There are examples given to show the im- 
provement over some of the other methods. 

Non-zero initial conditions are treated in Chapter V. It is 


shown how these become equivalent to an additional input function, and 
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can easily be incorporated into the Volterra series. Using these results 
the case of a step input is re-examined and a solution obtained by trans- 
lation of the equation which gives a better convergence than than given 


by the straight equation. There are examples worked out and a compari- 


son made. 
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CHAPTER II 


MATHEMATICAL PRELIMINARIES 


2.1 Introduction. 

In this chapter certain mathematical tools which are used in 
later chapters will be outlined. The object is not to give the rigorous 
mathematics behind these tools, but rather to introduce them in a 
descriptive fashion so that when they are used later there will be no 
loss of continuity. 

Firstly the algebra of working with operators is discussed. 
This is important in order to express a feedback system in terms of its 
constituent parts (plants), and is also very necessary to simplify 
the manipulation of terms of the Volterra series. Then single input 
feedback systems with one nonlinearity are considered, and a general 
system is used which can be reduced to several types of systems. The 
operator equation of this general system is the one that will be 
analysed in the remainder of the thesis. 

Norms play a large role in the functional analysis approach 
to the stability analysis of nonlinear systems, and the choice of the 
best norm to use for a particular system is an area of research that 
has not been looked at very deeply. Two norms are discussed and will be 
used later along with various methods of describing the size of the 
effect an operator will have on a function. This is called the gain of 


the operator but has also been called the norm of the operator. Several 
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examples which will be important are introduced. 
Finally the special class of operators which are contraction 
mappings are introduced and a fixed point theorem is stated along with 


some boundedness properties of this class. 


2.2 Operator Notation. 


A physical system which has inputs as functions of time and 
gives outputs as functions of time will be mathematically modelled as 
an operator. In this context an operator is a transformation or 
mapping from a space of functions of time into another such space. For 
this thesis only operators mapping a space into itself will be dealt 
with and the space in particular will be a normed linear complete space 
(a Banach space). 

Operators will be denoted by underlined capitals, such as H, 
while the functions will be denoted by lower case letters such as. x(t), 
y(t). Laplace transforms of operators and functions will be denoted 
by capitals which are not underlined such as H(s), X(s), Y(s). The 


fact that H operates on x to give y will be written as 
Y gic (e) = x) ix (251) 


As in the case of functions, the set of all functions x for 
which the operator is valid is the domain of H, while the set of all 
y which H(x) assumes is called the range of H. 

A numerical subscript to an operator denotes the order of 
the operator. For example, an operator of order n obeys the 


following rule where a is a scalar number. 


H (ox) = oH (x) (2.2) 
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A linear operator has an order 1. However, when the order of a 
system is mentioned this will refer to the order of the differential 


equation which describes (or can describe) the system. 


2.3 Operator Algebra [29 cy 

In order to handle interconnections of operators with relative 
ease, it is useful to look at the properties of the sum, product and 
cascade of operators. This helps in relating the behaviour of deters 
connections of subsystems (such as in feedback systems) to those of the 


components. 
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Fig. 2.1 Sum, Product and Cascade of Systems. 


The properties of the sum or product of two operators are 


obvious from Fig. 2.1. 
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(A BC = AG C) (24) 


The cascade of B following A is denoted by B*¥A. It is not a 


commutative operation but it is associative. 
(AXB)*C = A*(BXC) (2.55) 


The special cases of linear time invariant operators are 


commutative in the cascade operation. 





Fig. 2.2 Combination of Cascade with Sum and Product. 


The sum and product are distributive with the cascade as long 


as the sum or product follow the cascade as shown in Fig. 2.2. 
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The sum operation has a zero operator 0 which corresponds 


to the open circuit system whose output is zero whatever the input. 


O(x) = 0 
O+A=A 
A+ (-A) = 0 


(-A) is called the negative of A. 


The cascade operation has an identity operator I which 
corresponds to the short circuit system whose output always equals 


the input. 


I(x) 
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i*A 
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Provided it exists, the inverse of an operator when 


cascaded before or after the operator gives the identity operator. 


The operators used in this thesis can be expressed as a 


sum of operators of various orders. For example: 


where some of oie may be the zero operator 0. It will be 
convenient to consider that the zero operator can have any order. 


Consider now the effect of operations on the order of an operator. 
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Clearly addition has no effect on the order. Now if 


with ® a scalar number, then for 


Pp 

Similarly, 
Cp (ax) 

and so for cascades, 
Pp 


So it can be seen that 
sum the two orders while that of 


orders. 
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the effect of multiplication is to 


cascade is to multiply the two 
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2.4 Feedback Systems. 





Fig. 2.3 The General Feedback System. 


The feedback system of Fig. 2.3 must simultaneously satisfy 


the following equations. 


x4= iL s (e) 

gz, NG) 

yx L, (8) 

e=r-y (2.17) 


This is the general system which will be considered for the 
rest of the thesis containing a single nonlinearity N which is single valued 


and time invariant. L, and are linear plants either one of 


L 


which could be equal to kl a linear gain. 


Eliminating y, g and e from equations (2.17) gives 


x=L, (r-L,@@)) (2.18) 
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Now assume that x and r are related by an operator 
2 Silt 8) 
then substituting equation (2.19) into equation (2.18) gives 


Be) Se (ee Or) ))) 


Since this equation holds for all r it may be written in operator 


form as 


all 
Now suppose L. exists. Then, 


-1 * 
ZL, “H+ L*N*H = 1 


=2 


* * 
(L.* + L,*N)*H = 1 


-1 -l 
NX * * = 
(La + L,*N*L *L, )*H = 1 


-] 
+ L, *N*L )* *H = 
(1 LL N*L ) L. H =I 


-1 
= * ( + *kNX 


provided (1 + L,#N*L.) 7 exists. Compare equation (2.24) with the 
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(2.24) 


feedback equation of linear transform theory for the similar system. 


L,(s) 


H(s) = ; 
1+ L, (s)N(s)L, (s) 


assuming N to be a linear time invariant operator. 


(2.25) 
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So equation (2.21) is the operator equation of the general 
system under consideration. The output function could be any, Of @(t)\, 
x(t), g(t) or y(t) which can be obtained from equation (2.21) by using 
the necessary operators. For example if y(t) was the output and it is 


desired to express y(t) as 


Vee eG Gl) (2.26) 
then it follows that 

- = * * e 

G = 1, *NeH (2.27) 


2.5 Functional Representation. 


The solution and stability of linear time invariant systems 
can easily be determined by examination of the roots of the characteristic 
equation of the system. This method gives necessary and sufficient 
conditions for stability for such systems. In the case of nonlinear 
systems, the ideal would be to develop a theory which would enable the 
solution and stability of the system to be determined for any system. 

In this theory the linear systems would be a special limiting case. 
It is thought that using functional analysis techniques offers the 
best chance of achieving such a general theory. 

In engineering it is important to have some way of determining 
the size of errors or functions. This is taken care of mathematically 
by introducing norms [19]. The norm of x is written as ||x|| and mus t 


satisfy 3 properties. 


1. l|x]| antes |x| | = 0 iff x=0 almost everywhere. 
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2 I] x || = | a | heal where a is a real or complex number. 


3. [x+y I< Ixl+llyl , (2.28) 


There are no restrictions on how to choose a norm as long as 
it satisfies the 3 properties above, and so by judicious choice of a 


norm various properties of a function can be emphasized. 


For example, consider the problem of matching an analytic 
curve to a set of empirical points. If the error is considered to be 
the distance d. from the curve to the point i, and a measure of the 


closeness of the fit is to minimize max | d,| then there would usually 
L 


be a different matching of the curve than if () | ay ty? was minimized. 
af 


One of the more useful norms for functions used most 


frequently in engineering is the uniform norm [46]. 


|x || = sup] x(e)] (2.29) 
e 


The space C[a,b] of continuous functions on the interval [a,b] with 
the uniform norm becomes a normed linear complete space and so a 
Banach space. The space C[a,b] will be the one used throughout this 
thesis. When the uniform norm is used it will be used without sub- 
SGRIpE. 


Now consider the linear space C[a,b] and the norm (see [6]) 


xl], = sup e°" | x(t)| (2.30) 
te[a,b] 


This norm is equivalent to the uniform norm in that 
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Since the norms are equivalent, sequences converge in the A-norm 
if and only if they converge in the uniform norm. Thus the space C[a,b] 
with A-norm is complete and so is a Banach space, which will be called 
C, [a,b]. Whenever the }+norm is specifically used it will have the 


subscript As. 


2.6 Gains of Operators. 


A system which is bounded-input bounded-output stable has 
a bounded output for every finite input. If the system is represented 
by an operator H, then for the systema be bounded-input bounded- 
output stable implies that the gain of H is finite for all finite 


inputs. So if H maps functions x«C[a,b] into yeC[a,b] then 


Hx (2:32) 


< 
It 


II wl 


| Hx || (2.33) 


where for this discussion either norm can be considered. 


Consider the Lipschitz norm of an operator [19]. 


Definition: The maximum incremental gain of an operator H over a 
set S c D where D is the domain of the operator is defined to be 


the Lipschitz norm of the operator. 


inel| u]|= sup{ll¥Gc) - HG )I/ Ix, - x) ll 2 x, # xe Sh 
(2.34) 
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Definition: The slope of the operator H is defined to be 


Selim HG = Hx) / Go x) SOY 


dx X.+>+ X 


L 


provided the limit exists. 


Clearly it can be seen that for time-invariant operators 








inc|| H||= sup = 
x (2.36) 
dx 
Now if it is assumed that 
Oe =30 (2537) 
then setting x, = O in definition (2.34) gives 
JED = ound] ux I/Ix, Fs x, # 005) (2.38) 


and this will be called the maximum gain of the system. 
Since equation (2.38) was derived from equation (2.34) by 
fixing Xo» then the supremum in equation (2.38) is less than that in 


equation (2.34), as it is taken over a smaller range. Thus, 


inc] Hl] 2] HI (2.39) 
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From equation (2.38) it can be seen that 


Bx lls | Bll [ell ¥ xes (2.40) 


Substituting this in equation (2.33) gives 


ly < HI [Px || (241) 


and using equation (2.39) a coarser estimate is obtained 


lly || <  inc]H]| ]} x|| (2.42) 


2.7 Examples of Operators. 


(1) Linear Time-invariant Operators. 


For this broad class of operators the following applies 


ine | f:|| > = esup { LG) - L(x,) | / || x Xy || 3 x, # X2 eS} 
= sup { LG, ~ x») | / || x67 x | 3X) 7 Xy #0eS } 
(2243) 
inc ipa = || u|| (2.44) 

Consider the case of a linear convolution operator, which 

is the lined oneeainar to be used here 
y(t) = L(x) = of? n(tet) x(t) dr (2.45) 
with 

<Ct) 2 = htt) “=. 0, Fe<_0 (2.46) 


ei 


MmmA 
(Ve sa) 


(ia,S 


bend 


(EE.S) noftsups at alia gutivitsedsua. 


4 z tad 
BIAq trstirsaviatr—smil snout (f) 












} os 


; : ; ie. ii 
meas od mBo 3£ (BE. §) nobdeps - mor 


om 


WoW ite. to ~ eed 
Heel} BW > Way 5 
7 


ve 





Isetboo 6 (€f.8) nokisups gakau Bae 


erots1eq0 10 astqmexd 


we TE ce ne 








zwatlo beotd elf 257 


be 


: Nes mn 
a ee es 
SUSE ee h! 185 


20 


ly(t) | < (fi [n(e=e)| [x(r) Jar 
sup|y(t)| < sup|x(t)| sup, f° |h(r) | dt 
Ls t t 
yl Well [perder (2.47) 


It has been shown [15], that for these linear operators 
lz = (fo [nce) [ar (2.48) 


Thus this class of linear systems is stable for all operators which obey 
the following. 


af ln@)lar < « (2.49) 


Using the A-norm, the following is obtained from equation 


(2.45). 
ly(t) | < of |nct-1) | [x(r) |e" oo dr (2.50) 
=\T Eat 
< #sup ‘|x(r) Je sup |n(t-1)| [~ edt 
Te lO7>o] (t-t)e[0,T] : 
oe | 
= {x|,e cup gin (©) | — = (2.51) 
tealo sit 
ly(eyfe** < |x|, sup |nce)| + 


te[0,T] 


Then taking the supremum of the left hand side over te[0,T] gives, 


llvl, < Ix, sup [ace] 2 (2.52) 
te[0,T 


OF v : 
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(2) Nonlinear no-memory Time-invariant Operators. 


For these cases, the operators can be described by a graph, 
and the maximum incremental gain is the largest absolute value of the 


slope of the graph. For example, consider the nonlinearity 


P 
N(x) = ) n, (x)? (2.53) 
jal? 


Then using equation (2.35), it can be seen that 


iP 

dN i-] 

ee Gs (2.54) 
j=l? 

and so from equation (2.36) 
P aa 
inc ||N|| = sup | ) n.j(x)? | (2.55) 

xeS 'F=1 J 


From equation (2.38) 

















P a P hes 
|x] < xl} 2 n, (x)7 < ||x|| sup | ) n, (x) | 
jae xeS '4=1 J 
° j-1 
Issel / llxll s sup | Zn,o037* | (2.56) 
xeS ‘j=l : 


(3) Nonlinear Convolution Operators of Order n. 


These are of the form 


HG) = a een h (t-tyoeeestot )x(T))---x(t dt)... dt, 


(2557) 
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It can easily be seen from this equation that has order n. 


Consider now a bound on H (x). 
1 


|H (x) | 


| A 





tS & 
a | Pe Greer parent ele): P| 26 0c #) -+-[x(t)|dt,...dt, 


[A 


n co co 
xl of Bee Jb (tyseeot) [dt y.- dt, (2.58) 
Taking the supremum over t of each side gives 

a Gx) |< WPL FU tye eet ldty ee edt, (2:59) 
Thus for absolute convergence of the series given by 


H(x) = ) B(x) (2.60) 
j=l? 


it follows that 


IHG) <2 IH, Goll 
int | 


=.) IP Soe fds yee et) ety. edt, (2761) 


A sufficient condition that the series in equation (2.61) 


should converge, is that || x || be small enough, as long as 
age 2. fen (e)3- 22th) dey. - dt See = (2.62) 


The computation of the gains of these operators is more 
difficult than for the other types of operators considered, and it is 
better to try and find the gain of the total operator H by some other 


means. 
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2.8 Contraction Mapping 


Let E be a Banach space and let H be an operator mapping an 
open subset U of E into U; suppose there is a constant a , o<a<l, such 


that for every Xp »Xy say 1a) 


| H (x) - HG) < all x, II (2.63) 


or using the notation used here 


inc ||H|| = o (2.64) 


then H is called a contraction operator on U. 


If it can be shown that an operator is a contraction then 
there are many theorems available [20][21] to estimate the fixed point of 
the operator,and the error bound of the solution obtained by iteration,as 


well as a bound on the solution. 


Theorem [21] Let H be a contraction on U in a Banach space E with 
inc lle || = Hie Suppose there is a sphere S:{|| x-x|] < a} that is contained 
in U and that 
- < (1- 2.65 
I| Hx.) - xl < G-aa (2.65) 
Then there is a unique fixed point x~ in S, that is, a function x” for 


which 
H(x) <=; (2.66) 


and x~ is the limit of the sequence { x, } generated by 


x 


Sy LD (2.67) 
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and the estimate of error is given by 


[fran a x || < a™(1-a)~ "| x, - XII (2.68) 


Proof: (See [21] p.37). 


Comment: If the operator is a contraction over the whole space E then 
clearly a is as large as needed and so the condition given in equation 
(2.65) does not give any restriction. Also equation (2.65) can be 

satisfied by choosing a such that the left hand side is small enough. 
This condition makes sure that during the iteration, the iterands are 


always contained in U. 


If the operator H is in an equation such as 


x(t) = y(t) ot HG) (2.69) 


then it is possible to bound x provided H is a contraction. 


Hail < yi] +. [| B= 


by the triangle inequality of norms, 


[EAS Rod Sei ies) (2.70) 


by the use of equation (2.40). Then rearranging equation (2.70) gives 
-1 
H*ll <]] y] G-|[alp om) 


provided 


|HI| <2 (2.72) 
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But since H is a contraction it follows that 
pe feseeinc || ae [|= <1 (2.73) 


So for an equation such as (2.69) where the operator H has been 
shown to be a contraction in a particular norm and over a space then 
a bounded function y(t) implies a bounded solution x(t). This then is the 


basis for the bounded-input bounded-output type of stability. 
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CHAPTER III 


UNIQUENESS OF THE VOLTERRA SERIES 


See Lntroduction. 


The Volterra series is a part way solution to the problem 
of finding an explicit relationship between the output and input of 
a nonlinear system. It does give such a relationship but in the form 
of a (usually) infinite series whose terms are in general rather 
aifricule to evaluate Also the series is able to be formed only for 
a class of nonlinearities which can be expressed as a power series. 
There are two properties of the series which are important; does the 
series converge and is it bounded? 

Firstly, two methods of calculating the series are examined 
and it is proved that these two methods give the same series. Then 
methods are shown for evaluating the terms of the series and an example 
used to illustrate the computations and also to introduce the next chapters 
which will deal with the convergence and boundedness problems for the 


general class of systems under consideration. 


3.2 The Volterra Series by Iteration. 


Consider the system equation (2.21) with 


=. (3.1) 
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Then 


ay 


(3.2) 


This equation can be solved by iteration, assuming for the moment that 


the series converges and is bounded. 


Let 


be the starting approximation. Then 


HH’ = L_— LAN4&L 
= ae 
ig = L —- L¥N*(L_ - L&NAL ) 
= ye ee 


= ="LSN* (l= LAN*(L — LINX 
HO = L. - LAN*(L- L*N*(L.- LAN*L )) 
ol ie L - Laney” 
- Bab ee 


where H™ is the nth iterate. Now let 


M ; 
NX = i n(x)" 


j=2 


then substituting in the starting approximation gives 
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g=2 j=4 
where the H,* are given by 
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Substituting equations (3.7) into equation (3.4) 
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(3.8) 


M , 
1 
Now ( ) bo can be written in ascending order of operators, 


i=1 


where Le denotes the operator of order j[16]- 


Yd 
( ) Hs a= 
‘=r 
where yt GG) 
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yi > 
"aaa B 
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H, 


1(j 
i 


=] 


rf 


Mj 
ie 


i=Mt+j 


ite 
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P._o=t 
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Po 


gl cdae Panky 


is of order i+ j - 1 and is given by 


Love 


+1'~i-p 


(an) 


+L 


(3.10) 


1 A 
This expansion gives all the combinations of H. "s which when 


multiplied together give the Hi of order it+j-l. The os of equation 


(3.9) are operators of higher order of the expansion which cannot be 


expressed in the form of equation (3.10). 


into equation (3.8) gives 
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Substituting equation (3.9) 
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2 dened Leek L. 3 
H = Se. [ee H. H. +H. 4H te ay ( 


_ ame 2 3 

= 2n,L*(L, G* CL) ))-n,L*(L_) (37515) 
at Lo 4 11k) . 
Le ae a 2<j<M (3.16) 


There are Me of the a terms and by the nth iteration there 
would be M” terms. Some of these terms cease to change after a certain 
number of iterations, and so just add to the work without providing any 
additional information. For example, in the case above the terms Hy 
and Hy are equal as are Hy and Hy but the third order terms differ as 
shown in equation (3.15) where another factor has been introduced over 
that for Hy. 
3.3 The Volterra Series by Substitution. 


Consider the same equation (3.2) with N as defined in 


equation (3.5). Then 


M 
1 
H = L ead) n, () G17) 


G i=2 


Now let H be a Volterra series of the form 


Hee ) of, (3.18) 
jel? 


where Le is the nonlinear convolution operator of order j as defined 
in Section 2.7. Then substituting equation (3.18) into equation 


(3.17) gives 
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Now using equation (3.10) to express the sum raised to a power 
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(3.20) 


where bse has order itj-l. Using equation (3.20) to simplify the 


notation in equation (3.19), gives 


fore) M co 
af 
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a=) Wen pen? 
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HT] 
ie 
I 
i c-18 


by rearranging the terms and changing the order of summation. 


(3-21) 


(37.22) 


Now using 


the rules for cascade sum and product as outlined in Section 2.3 and 


equating orders on each side of equation (3.22). 


HoT 
min (M, j) ( 
k) ; 

= = * 
H, Lo REE} kt ea 
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= 2n 2 L4(L_ (L*¥(L )2))-n L*(L je 
3 Nae: eg es a ae 


By comparing equations (3.25) and (3.15) it is seen that they 


are equal. However will they remain equal through more iterations? This 


is answered in the next section. 
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Evaulation of the terms of equation (3.23) give as many terms 


of the Volterra series as is desired. 


3.4 The Uniqueness of the Volterra Series [23] 


Consider now what has been done in the previous two sections. 


For the iteration case an initial estimate was chosen for the solution. 


Let 

oe Cr) 
and then 

x = L(r) - LAN(x) 
Now take the initial estimate to be 

0 ¥ L(x) 


Then from equation (3.27) using iteration 


x, = Lr) - LAN) 


- L* 
saeiay L,{) LAN (x) 


So that the limit of the sequence {x } 


= 1 = ma | * 
x lim x L, (4) lim LAN (x) 
m0 no 


For the substitution case it was assumed that 


n 
x” = lim ) H, (r) 
mo j=l 


assuming the limit exists. Substituting in equation (3.27) 


n 
f Lit) LAN Ca 14 Hk) 
ers lin 2% 


L(t) - LAN(x’) 


(3.26) 


(27) 


(3.28) 


(3.29) 


(3.30) 


(3532) 


(3.32) 


(3:33) 


(3.34) 
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So it is necessary to show there exists a unique solution to 
equation (3.27) so that 


lim L*N(x ) = L&N(lim x ) (3.35) 
no«o uy n>co R 


Now if L*N obeys a Lipschitz condition that is 


I] LANG) - LANG 4) || < Rll x, - x (3.36) 


Sail 

where K is a real positive constant, and if K<l over a region U, then the 

contraction mapping theorem of Section 2.8 can be applied to the system. 
So there exists a unique solution to equation (3.17) if the 


conditions of the theorem hold, and consequently equation (3.35) is true. 


It has been shown that the first few terms of the two series 
calculated by the two methods are in fact identical and it is found that 
as the iteration is carried on further,more terms do equal those calculated 
by substitution. The problem with the iteration method is the large number 
of terms which must be carried through each iteration. The easiest method 
of evaluating the series is to use the substitution method to get terms 


up to a certain order and then perhaps do one iteration. 


3.5 Evaluation of the Kernels of the Volterra Series 
et Se Aerneis or the Volterra series 


There are two main methods of evaluating the kernels of the 
series. The first uses time domain analysis by using the integral 
expression for the general convolution operator. Then the evaluation 
of the kernels involves integrations. For example consider equations 


(3.24), (3.25) with 


LL) = ei ECCT rat) da (3237) 
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L(r) apa e(t-t) r(t)dt (3.38) 


and 


te) 


2 (tae) ome Oe mer orart <0) (3.39) 


as is found in all practical systems. Then it can be seen that 


H(t) = ef F(t-t) r(t)dt = NCES phe ae fale, 


(3.40) 


HG@)-= of f° hy (ome tat hor (x, E(t) drat, 
(3.41) 
= ony (f° g(t-r)dt ((f* £(t-r,) x (r,)at,)° 
| (3.42) 


= “n, Sf ese g(t-T) Gao) E(taas EG) 


r(t,)dt,dt,dt 


Now using equation (3.39) the integrals can be expanded 


Hye) = -n, [TS fh alter) £(e-t)£(r-t)dt (t,x (ty) dt dt, (3.43) 


Thus it can be seen that 


hj(t-t,, t= = -n, of B(t-t)£(t-1,)£ (t-1,)dt (3.44) 


>) 
and similar integrations can be used to evaluate the higher order 


kernels. However, as in the linear case, the integrations become some- 


what easier if they are done in a transform domain. 


.. 
EE ; 
_ ° 
~ 9 
P ) 
(8é.£) » sbCHY Gade” | = Gg 
bas 
_ 
(ee .€) O>4 voi O= (+ = ee GR J ; 
Jed mesa sd aso Jt med? yehedaye Laoticarq Cie ‘at bawe® et es 
rb(,7)2¢,t=2)d°° | = rbCpe (re3)a |. * Ge E. 
pg teehee eS: eae - 2 : 
(Oa > 2B, : . 7 2 
h thé ( ~ oh ie a ).H ; 
gM gOS AG 7) t Corn ts TEE hae oe ae 
([8,€) g 
Y b(.7) x =3)% | )- ThAT=a)g 7 i = 
‘f Gg : ae ee” Se 2. J 
($8) am 
( rir = t( TT)? (r-3)2 ‘4 Ty 4 a 
ed i az Lo “Oo *o $$ 
ae 
| : 
rb, Tb, THC. t)4 | 
<° =a 
hbsbaayxs $d nso alssgotat sds (C8,£) paotzaups anias wow a 
oy 
4 id 
(€2,£) <shestites) as aioe Tae he = Gel 5 





hecbedembestes: 
edi Sie ule i a 


ss (739 , see 


~~! 


34 


The Laplace transform in the linear case changes convolution 
integrals into algebraic multiplications. For the nonlinear case, 
multidimensional transforms have been introduced [11], which allow the 
nth order convolution to become an algebraic multiplication. These 
multidimensional transforms are defined in a similar way to the Laplace 
transform [8]. If F(t,,---.t) is an nth order kernel or function and 


F(sy,-++s8)) is its multidimensional transform, then 


co oO 


= Gm ts wero SHne 
F(S)2S59++++58)) = Ted F(t,,-+-,t Je ‘Loot n Ndt,..dt_ 
(3.45) 
Beware (ee ea) tt 
J.” n 2045 , a: an 
o,-j2 oa -jx 
eemS Li 

ds,--+-ds | (3. 46) 


Many of the properties of the Laplace transform have similar 
properties in the multidimensional transforms. A few of the more useful 


properties to the present application will be listed. 


Convolution: In a linear system, convolution in the time 
domain corresponds to multiplication in the transform domain. This also 
holds in the case of a multidimensional convolution. Such a convolution 


would be given by 
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For this case it can easily be seen by taking transforms that 


NE gee Ue TG Rae Cpe cod CO (3. 48) 
This is not quite what is required as x should really be a function only 
of t but this can be accomplished in equation (3.47) by letting 
ChB et eo Et (3.49) 


This is known as association of variables, and for rational 
transform functions it is possible to perform the operation in the transform 


domain by inspection, with the help of a table of examples [11][25]. 


Cascades of Systems: Two types of cascades will be used: 
a) A linear system L, (s) followed by an n-dimensional system with 


transform H (syo+++8))- The transform of the cascade is 
* = 
Hq, Ly aa H (Syoee++9S,)by (S))+++ Ly (s)) (3.50) 


b) An n-dimensional system H_(s,>++++98)) followed by a linear system 


L,(s). 
L*H — L,(s) oy +. ..ts) H(Sy2+++2S)) (3.51) 
Now consider equations (3.23) again 
Hr 
u, ite isin ee (3.52) 


Let the transforms of L. and L be Ls) and L(s) respectively. Let the 


transform of 4 be LAH Sod taHe Then it follows from the above that 


Hy, (eyee L,(s) (3.53) 
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H, (sj 585) = ‘st, L(sj+s,) L(s)) L_(s,) (3.54) 
H3(s)>S5s83) = -2n, L(s,ts,ts,)H, (s,)H,(s,.84) ~ 


nl (s)ts,ts,)H, (s,)H, (s,)H, (s4) 


(855) 


s 2 
= L(s,ts,ts,)(2n, L 


a (5, )L(syte,)L, (s))L, (84)- 


nL, (s,)L, (s,)L, (s4)) 


H,(s)>S5583) = L(s,ts5ts,)L, (s,)L, (sq)L, (s4) (2n,L(s,+8,)-n4) 


(3.56) 


The expressions become quite lengthy as the order of the kernel 
increases although the basic simplicity remains. However, the association 
of variables where required, although only needing simple steps, becomes 
quite cumbersome. 

Once the transforms of the kernels are calculated then many 
types of inputs can be treated including classes of random inputs [25]. 
This method is a generalization of the method of Poincare(see [12] p.198) 
and also that used by Doetsch [8], in which the solution is found for a 
particular input. For sinusoidal inputs the transforms of the kernels 
can be used to give the steady state sinusoidal output of the system and 
the harmonics [11]. This can be useful for distortion analysis of systems. 
Obtaining the response to a step input of some magnitude can give the 
speed of response of the system and the series is quite accurate for the 


first transient response of the system. 
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3.6 Example of Volterra Series Calculation 


Consider the general system of Fig. 2.3 with zero initial 


conditions and 


i ae (3. 57) 
a 3 
N(x) = e(x) (3.58) 
L(x) = L(x) = ft oR xcryat (3.59) 
Then the operator equation of the system becomes 
3 
aCe) e (ee felis) (3.60) 


Now using equations (3.57), (3.58), (3.59)with equations (3.23) gives 
= els, 
j22 (3.61) 


Now assuming that 


pe =o 
=o — 


and grouping like terms together for convenience it follows that 


= Ba6Z 
H = 0 (3.62) 
H, = - au) =— cua@)>=- eG)” (3.63) 
= 3.64 
a, = 0 (3.64) 
. eye 2 
H. i ea eL*H, so EL*3 G,) Hy 


z (a: veneg i) > TSE CL)” 


= 367L* [(@)7*@)*)] (3.65) 
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Expressing these operators in their transforms gives 


| 1 
H, (s) = aa (3.68) 
Hy (Sy202 99 Son)= 0 (3.69) 
Hee (se 3S... 1S.) = -€ a a (3.70) 
Gere] 2-33 s,ts,ts,ta (s, ta) (s,ta) (sta) 
2 
3e€ l 


He(s..,S), 8.46758. = { Pa wal es. A. 1 ae Le 
Bee eet a (s,ts,ts,ts,ts ta) (s, ta) (sta) (s,ts,ts +a) 
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rb) se ke u(t) (3.73) 
R(s) = = (3.74) 
R (s). 5s.)= oH (s)5+++5S,)R(s,)----R(s)) (34.75) 


Now using equations (3.55) to (3.59) 


So 


2, (s) ica) (3.76) 


Eko le 
s(sta) _ 


t 


Z(t) = © (u(t)-e 7°) (3577) 


This is the linear approximation to the system. 


a 























Xis.s,s,) = -£ ®— eae -jape SS 
Sie he Be ee + 
a ‘a (s) s,ts,ta) Sy s,ta Sy s,ta 
tee sta 
ody CE aa ORE ike 
aa s,ts,ts,ta S48 s,(s,ta) 
1 BE ft } 
Sere ee eee eer | ee ).(3.48) 
s,(s,ta) (s, ta) (sta) 8, sta 
Now using the inspection technique to associate variables, 
first Sy and S5 are associated. 
h, (apes. 5. il = - £534 ¢ - 2,15y¢-4 
J glee ee s,@s, a‘a s,ts,ta Sy) s,ta s,t2a S3 s,ta 


(3.79) 
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Associating Sy and S, now gives 




















£.(s, 484) --£*345)2¢-34,—4._-) 
Be bee | Sy @®s, a ‘a’ ‘s,ta’ ’s sta s,t2a s,t3a 
(3.80) 
A 5 POTS AY ha Serie eee te ee 1/2 
ae (s) a = [ - sta (sta) 2 st+2a ue ar38) 
(3.81) 
A k. 3 a =. es 
xX, (t) = - ea) (u(t) te at - 3at)-3e “+e ao 
(3582) 
In a similar fashion, 
Q.(s s.) = 35,24 ,—4—_, iT - : 
St ela aa 5 s.ts,ts_ta’.'!‘s. sta 
3° \4e05 i=l “i i 
() s,)+a 
, i 
i=1 
(3. 83) 
and associating variables gives 
A ee Vex dyik Sia en ees 3) HE 
X,(s) = 3©) oe sta Ss aa sta (sta) ° st2a e (sta) 7 
8 3a 4 igg: 
st3a  (st+3a)? st4a  st5a } (3. 84) 
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(£) Actual Solution 
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2(d) (a) Linear approx. xy 
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(b) xy + X, 
A A A 
(c) x + X. + Xe 
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Curves No. 13 G] = 0.5 
(£) Actual Solution 
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(a) Linear approx. xy 
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(c) x + X, + Xe 
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Time t (secs.) 


Fig 3.1 Approximations to System 
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and similarly for x except that here the expression must be split into 
a sum of two simpler parts in order that the association of variables 


may be carried out by inspection. 


%y(t) = ~ (2)? "(120 (t)+e°#* 26. 737546. 375at-11.25 (at)? 


-3at 


Pate S4atas4 (at) ice 


Aes Gaeyal ce [36. 5625+ 


78.75at-20.25(at)“]-e 19*[23-48atl+e 9** [22.3125 


6 


SGey3 1556 


25. 625a then see alat) 


(3.86) 


These results are plotted in Fig. 3.1 for the following parameters 


along with the exact solution 


a = ©) =1 (3.87) 


& 0.5 and 2 (3.88) 

It is seen that there is good convergence for small time 
intervals and the larger number of terms used the better the convergence. 
However for the Ee = 2) curve the solutions diverge wildly past a certain 
time while for the other curve convergence appears to be for all time. 

If the derivatives with respect to time are removed from the system 
differential equation, a nonlinear algebraic equation is obtained. Attempting 
to solve this algebraic equation by an iterative process such as Newton's 
method, will fail for larger Ss (such as Ke 2), and the successive 
iterates will diverge in much the same way as the solutions here have. This 


explains why the convergence at steady state is so poor in these types of 


examples. 
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CHAPTER IV 


A LARGER REGION OF CONTRACTION 


4.1 Introduction 


One of the major problems in the analysis of nonlinear systems 
is to obtain realistic bounds on the nonlinearity to insure that the 
system will be stable. Or the other way of tackling the problem is to 
find or modify the linear part of the system to give a stable system 
with the given nonlinearity [30], [31]. Each of these approaches has 


been used and the one to be investigated here is the former. 


A method is presented which gives good results for first and 
second order systems (meaning systems whose linear plant can be 
represented by first or second order differential equations), but does 
not appear to work for higher order systems. For these systems another 
norm is introduced which then gives the required contraction but does 
not allow any conclusions to be made about the stability of the system. 
However it does put a bound on the output for a specified time interval 


and so implies convergence of the Volterra series over that time interval. 


4,2 First Order Systems 


Consider the system of equation (3.2) where 


L(s) = rare (4.1) 
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<8 for all x « C[0,~) (4.2) 


where dand 8 are real constant to be determined such that the sector 
[a,8] is the largest possible for stability of the system. The system 


equation is 


x(t) = L(t) - LAN(x) | (4.3) 


*L =L | (4.4) 


In order to keep the system and analysis as general as possible 


let 


L,{) w(t) (4.5) 


x(t) w(t) - L*N(x) (4.6) 


Now expressing equation (4.6) in its differential form using equation 


(4.1) gives 


xk + ax = wtaw - bN(x) (4.7) 


Add bux to each side of equation (4.7) where u is a real number. 


x + (atbu)x = wt aw + b(ux - N(x)) (4.8) 


Now returning to the integral operator form of equation (4.8), 


x(t) = H(x) = w(t) - uK(w) + K*(ux-N(x) ) (4.9) 


where 


as 
K(s) = area es (4.10) 


) 
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The Lipschitz constant of H as defined in equation (4.9), is 
obtained as follows. 


N(x, )-NCx,) 


| HGx, )-HGx,) || = Re (a Geax) Gc =x) X17x 


2 
for all x, #x,€C[0,~) (4.11) 


N(x, )-N(x,) 





[HiGe) =H Gee) eK es x7 | sup u- 
A é in x x,€C[0,~) Xyeay 
(4.12) 
|K|| can easily be evaluated. 
i 20 ~(atbu)t ee b 
inet 2 dt = = (4.13) 


From equations (4.12) and (4.13), there exists a contraction 
if, 
: N(x, )-N(x,) 


ex sup a=) i 1 (4.14) 


TK 
x, #x,<C[0,~) elie? 


So, by removing the absolute value signs, equation (4.14) becomes, 


N(x, )-N(x,) 


a a 
is meer nel a eae < yu HF : x,#x,€CL0, ) 
Ab tye 
N(x, )-N(x,) 
wy + u(i-y) < —S «+ ul) (4.15) 
b = K47X5 oD 


If N(x) is assumed to be differentiable, and letting y>1l 
with the consequential removal of the equality signs, equation (4.15) 


becomes, 





< 2Qut— (4.16) 


Since u can be made as large as is needed, the following condition is 


obtained. 
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dN 
= SE ¥ xeC[0,«) (4.17) 





For this system, these limits [- »©) define the Hurwitz sector,over 
which the system with N(x) being a linear gain, is asymptotically stable. 
Since the linear case is to be included within the general nonlinearity, 


then this sector is the maximum possible for stability of the system. 


Now consider a bound on the output of the nonlinear system. 
x(t) = (1 — uk)w(t) + K*(ux - N(x)) (4.18) 


Taking norms on both sides and rearranging gives, 


N(x) 
ex's 


lel} < [[a-uxll [wll + (Kl sup IIx|] (4.19) 


x#0eC[0,~) 


u- 








It has been shown that for N obeying equation (4.17), and 
for sufficiently large u, then 


N(x, )-N(x,) 


IKI | sup deere ee ee (4.20) 
x, #x5€C[0,~) bs 92 


and along with the assumption that 











N(O) = 0 (4.21) 
these imply that 
Ki heup) Nite tee ee alli a sup eee sare he 
x, #0<C[0,~) i x, #x,€C10,~) ene 


Y Cx al (4522) 
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Substituting equation (4.22) into equation (4.19) 











lx |< I] z- ok] Iw] + y [xl 
lx I< Gey" a - xl Iw I (4.23) 
<  Gey)* I 2- oki ll I (4.24) 
|| Z - uk]| <1 + te ee ae =1+ a 
SZ | (4°25) 
Jol < max (itaeat seh ee (4.26) 


by the definition of the system. 


Clearly then equation (4.25) implies bounded-input bounded-output 
stability for all first order systems with a single nonlinearity which 


obeys 


vexe eC. [0.2) (4527) 


4.3 Second Order Systems 


Proposition: Consider the systems which can be described by the 
operator equation 
x = L,@) - L*N(x) (4.28) 


where the Laplace transform of L is given by 
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L(s) = 2 ee (4.29) 


and L, is a linear plant of order less than or equal to 2 such that 


a <2 (4.30) 


and nonlinearity N(x) is defined by 


M A 
N(x) = 2 n, (x)J (4ca1) 
1 
dN M + 
a a = (is in, a)” : = 6 (4.32) 
4=2 


where [a,8] is the Hurwitz sector for the linear system with N(x) 
replaced by a linear gain; then equation (4.28) has a unique solution 


and is bounded-input bounded-output stable. 


The Hurwitz sector for second order systems, with all the 
coefficients in equation (4.29) positive, is the sector [a,~°), where 


a, a 
a 2 ming Cae. ea) (4.33) 
2 


If some of the coefficients in equation (4.29) are negative, 
then the 8 of equation (4.32) may be constrained to be less 
than infinity. 
In the proof of the proposition, complex valued functions will 
be introduced. Their absolute value will be taken to be the square 


root of the sum of the squares of the real and imaginary parts of 
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the function. 


EROOE: Define w(t) to be 

w(t) = L(Y) (4.34) 
Then equation (4.28) becomes 

x(t) = w(t) - L*N(x) . (4.35) 


and expressing this equation in its differential form gives, 


age e — ° ids 
xta, xta x = wa wta,w - (b, aS +b, N(x) (4.36) 


Now add to both sides (2e-a, )#(2e7-a,)x; where c is a real positive 
number. 
forsh2e x = wra.wta.w + Ee Oca )x-b,N(x)] + Oe a ) 
be deal dt if I; Z 


x - bwNCx) (4.37) 


Define two new variables. 


Ky = x, + (2c-a,)w - [(2e~a,)x-b N(x) ] (4.39) 


then the following matrix equation can be written. 


(4.40) 
The system matrix can be diagonalized using the Van der 


Monde matrix. 
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The eigenvalues of the matrix are -c(1+j), where j=v-1. 





' 1 1 
- ee Pash | Coal 
1 . 
= Ed) 
Bee tak i (4.42) 
g(lt3) 2c 
Now transform equation (4.40), using 
x 2 
& = (| (4.43) 
2 2 
Z4 -c(1-j) 0) (24 $(b,-5([b,/c]-b,)) 
Bll Ile) - [3 Bi i 
2 0 -c(1t+j)? “2 5(b, +4 ([b,/c]-b,)) 
$(2e-a,-j(a,-Lay/e])) | 
| 1 (x-w) (4.44) 
5 (2e-a, +5 (a,-[a,/c])) 
and from equations (4.34), (4.41), and (4.38), 
x ite SE z,+z. (4.45) 
dN (x) 


Substituting equation (4.45) into the expression for a5 





given in equation (4.32), 


dN (x) M + 
rl a = ) jn, (2,42540)” : (4.46) 
a2 
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It can easily be seen that 


dN (x) a dN (x) dN (x) 


dz, Zo dz., Z4 dx 


const. const. 








(4.47) 


Consider equation (4.44) broken into its two separate 


equations, and consider the equation in zy first. 


Z4 = 5(2c-a,-§(a,-[a,/e])) (2,-w)+ $(-a,-5 (a,-2e-La,/e])) 24- 


$(b,-1([b,/e]-b,)) NC) (4.48) 


Now define real constants i A 


> de> Myo and add (e-jf)z, to each side 


of equation (4.48), where e and f are real numbers, such that 
Z +(e-jf)z = t(2c-a -j(a,-[a /c]))(z -w) + (e-A -j(f£-A jz - 
ay a Gi: 2 1 1 2 2 1 2 ih 
O75 NG) (4.49) 


and putting equation (4.49) into operator notation gives, 


2, = H(z,) = $(2c-a,-j(a,-[a,/e]))K(zy-w) + K(Le-d,-J (£-A,) J 


24-(4-5d,N(x)} (4.50) 


Finding the Lipschitz constant of the operator H gives, 











[iée, =H) Il 0e[ RAP [leas =e. ee sup eae=") Cs 
ee it cae ae ie z,€C*[0,~) > ct) 
5 (f-r,- A oe) (4.51) 
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where C*[0,“) is the space of continuous complex functions. Then using 


equation (4.47) and noting that aed) 


ay 


is a real function as well as 


the following 


| K |] = af lexp (-etjf)t |dt = of exP (-et)dt =< (4.52) 


gives 





dN (x) dN(x) .2 
3 y+ (EA, 1, BO 22 


ll 243 - #22 (4953) 


For contraction to hold for H, it is required that 














al we ee dN (xe? i AN) eas Z 
- a Coty prem e a Gants Viera @ JZ < <1) (4.54) 
and this implies that 
dN 2 2 2 dN 
LOD)? Py 29-2) yleapreng ag) Pee (ey) EAD)” 
for all Z4 —eC*[0,~) (4.55) 


Completing the square and taking the square root and rearranging gives 
Dy eet peer 2 2 2\r 


a 





2 a Zz 2 
<Og +, @: <),(e-¥,) +4,(E-%)) + es ee is 


(£-,)°) +O (eA) 42, CE 7 F (4.56) 
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Letting +l and removing equalities makes equation (4.56) of the form 


dN (x) 
ae < B (4.57) 





Now if N(x) were a linear gain, then for stability of the 


linear system so formed, it would be required that 


a, a 
a=- min ————) (4.58) 
1 


and substituting this in equation (4.56) gives 


2 2 2 
7 f4cb,t+(e-f) 4c b,+2(b,-b,c) (a, c-ay-2e )-2a,b,¢ 
oh dd ae 2 
2b) c+ bo - 2b,b,¢ 
(4.59) 
b a 
m ay ae ee OR al 2 

e(a,+ ab, ) = ( f-ct+ 9 (a, tab, on =») +- i (a b,+a,) | (4.60) 


The equation in Zo is just the complex conjugate of that in Z4 
and so can be treated in like manner to give the same conditions; namely 
that contraction holds for all N(x) obeying condition (4.57) with a and 8 
defined in equations (4.58), (4.59) and (4.60). So by the contraction 
mapping theorem there exist unique functions z4% and Zo» satisfying 
equation (4.44) and so there exists a unique x(t) defined by equation 


(4.45) which is a solution of equation (4.28). 


Define 


Niza) 8=) (Bests (f-A,)) 2, - Cy-5,) NG) (4.61) 
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2 SE ES (CP eC) . (4. 62) 
In a similar way 
2 oyete zon Vo ue (z,-w) +K' AN" (2.5) (4.63) 
where v, and v, are complex numbers. From the previous results on 
the contraction it is known that 
|| KAN' |] < ine || KAN" || < y,< 1 (4.64) 
| K'#n" |] < ine || Kt#N"|] < vy, <1 (4.65) 
Taking norms in equation (4.62) 
2, lls Ivy) WS ew + Pk ell ay (4.66) 
lz, Il <ly, [ell Geyp 2c! + Iw I (4.67) 
Similarly from equation (4. 63) 
zy Ws Ivyl I'l G-y.)* Clle,ll + Well (4.68) 
Now substituting equation (4.68) into equation (4.67) 
lz I< lv, tx Gey) wll a+ |yy} Ix @-) + 


-1 -1 
Iv, volIKINK' IH G-y)> Gey) ~ Thy (4. 
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Using equation (4.52) and also 
|K'|| = 4 lexp(-et-jft) | dt = of exp(-et)dt + (4.70) 


gives 


Iz, 1 (1-|v,v,JeEd-v,)G-v,)1) < Pilwl| (arn) 


where D is a finite positive number. Now e can be made as large as 


necessary, so that there exists an E, such that e2E implies; 
| 1-|vyvp|[C-7,) Q-y,)17* &? > 8 > 0 (4:72) 
Thus equation (4.71) and equation (4.72) give 
lz, ll < 6D Ill] (4.73) 
and putting equation (4.73) into equation (4.68) gives 


llz,ll< Ivy[feQ-v,)1> +8) Hwl] 7a) 


Equations (4.73) and (4.74) imply bounded-input bounded-out put 
stability for the system described in the proposition. It can be seen 
from equations (4.59) and (4.60) that the sector [a,8] can be expanded 


to the Hurwitz sector by increasing e sufficiently. 


Examples of Sector Calculations 
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and this implies 








= 2 Cy Ra eee (4.75) 
by dx b, b, : 
<I +2 
and clearly as e+», so does (2— f - <>) 
b b 
2 2 
a 
us 
Bye b. = 0% seo ae = — 
2 
eh 
f 52 
Let c =¥ -—> then 
2 
f =v 2a. 
and this implies 
EP “gdN(2)? > a i 
< < = i 
a dx Bie wane aes 


and this, as before, can be expanded to cover the Hurwitz sector. 


4.4 Higher Order Systems [24] 

The diagonalization procedure used for the second order systems, 
and which was a generalization of the method used for the first order 
systems, does not appear to give the larger region for third and higher 
order systems. For these cases it is necessary to introduce a new norm, 
the \- norm, and the space C, [a,b] of continuous functions on [a,b] will 


become the space Cc, [0,T]. 
Consider the same system equation as before 
x = L,@&) - L*N(x) (4.77) 


with 
aU < © (4.78) 
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L(s) = rr m>3 (4.79) 

s ta_s ata +a 
al m 

s ik 

N(x) = n, (x) (4.80) 
i=2 

5 dN(x) _ oo Atal 
a rea BL noc? =.8 (4.81) 
x(t) = w(t) - L*N(x) (4.82) 


Add to each side of equation (4.82) (c-a_)L(x), where c is a 


real positive number. 


(It+(c-a_)L)x(t) = w(t) + L((e-a_)x-N(x)) (4.83) 


x(t) = (I#(c-a )L) w(t) + (I+(c-a JL)" #L((e-a_) x-N(x)) 


(4.84) 
Now taking Laplace transforms of (I+(c-a_)L)"* gives, 
oy cTtas™ te. ta_ 
Cis Cota Es Ce.) = a (4.85) 
= s"ta,s” “+. ..+c 
= (I#(c-a_ JL) *(s)L(s) = ———+ = K(s) (4.86) 
s +a,s aro Goare 
x(t) = H(x) = (I#(c-a_ JL)" w(t) + K((c-a)x-N(x)) (4.87) 


The Lipschitz constant in the \-norm is now to be determined. 


|H(x,)-HG,)| = [K((c-a,) @y-x,) - ING y)-NG,) 1) (4.88) 


where K is an integral operator with kernel k(t), as given in equation 


(4.86). 
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|Hx,) - Hx) | 


JA 


af leCt=t) | | Coma) Ge.) (1) = ING, )-NGx) J] dt 


(4.89) 


A 


sup Jkce) | J sup | (esas) Go ox)) - [N(x )- 





























te [0,T] tel, i 
N(x,) ] [dt (4.90) 
~ dN 
sik f sup |x 17%, An eee sup e-a = a 
te(0,T] xeC, [0,7] 
dN (x) At_y 
< K ||x,-x,]| sup Coa eva, < 

tN kN 4 C, [0,7] m dx r 
(4.91) 

where K is defined to be’ sup LRGE) From equation (4.91), 

tefORT) 
rearranging gives the following. 
i - dN (x) 
|H(x,) - H(x,)le me oe : sup (hedtel a aey | x,-x, || 
= ee ene 7 m dx les uN 
xeC, [0,T] 
a dN (x 

—> | H(x ) - H(x vik < KA ‘i sup c-a _ ey Ix -x II, (4.92) 

=r: a2 —~ m dx Bead 

xeC, [0,T] 


For contraction to hold, the following condition must be 


satisfied. 





dN (x) 
Per dx 








< Y < 1 (4.93) 
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cca + (4.94) 


— C= =e. | 





Let y—1 and remove the equality and set = cK gives 


dN (x) 


th dx 





< 2c-a_ (4.95) 


Now for the inverse operator (I+(c-a_)L)* to exist, such 


that its norm is finite, c is restricted to a range of values such that 


m-L 


the zeros of the polynomial s'ta,s 


+...+c are all in the left half of the 
s plane. It is found that for equation (4.95) to cover the Hurwitz sector, 


c does in fact lie in the range to give a bounded inverse operator. 


Consider now a bound on the output 
25) : 
x ll, sll @+ (ea) bo] iw ll + ancl ull, Hl x I, (4.96) 
It has been shown above that for N satisfying equation (4.95) then 


ine [HI], <2 (4.97) 


> xl< ant }a+eay ol, ievll (4.98) 


Equation (4.98) implies bounded-input bounded-output stability 
in the space c,[0,T]. However this does not mean that |x(t) | is not an 
increasing function of time but that it is increasing no faster than 


De AE where D is a constant. 


By the contraction mapping theorem there exists a unique solution 
to the equation (4.77) in the space c,[0,T]. This solution can be obtained 


by the Volterra series which converges on the interval [0,T]. 
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4.5 Comments and Comparisons 


All the previous results hold if the nonlinearity contains a 
linear term. Consider the nonlinearity to be of the form nx + N(x) 


where N(x) is as defined previously. Then equation (4.3) becomes 


x(t) = L(x) - nb) - LANG) (4.99) 
(I + nL)x = L(x) - LANG) 


x(t) = (L+n,L)*L, (x) - + nL)” "LAN (x) (4.100) 


: <1" ds 
provided (+ nL) exists. It can be seen that (4.100) is of the same 


form as equation (4.99) with L, replaced by (1 + n,L)*L, and L replaced 


1 ; ' 
L) *L. As defined here, this transformation does not increase 


by (CL + n, 


the order of the system. 


Although the method used for higher order systems does not give 
stability of the output of the system, the systems do appear to be bounded- 
input bounded-cutput stable for the conditions imposed on the nonlinearity. 
This has been noticed by computing solutions to various systems. Also 
systems with zeros could be handled by this method but as stability bounds 


are not obtained these systems are not shown here. 


For comparison with other results, Barrett [2] considera the 


system 


x + a,x + ax + erty (r) (4.101) 


and found a bound on x for which the solution was bounded-input bounded- 


output stable. 


1 
|x|] < Gen) 2 (4.102) 
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where 
il: 2.- i 2.- i 
aasition [-a} (4a,-a, ) 2))/(l-exp [-a, (4a,-a] yew) 
H = 3 a}%< 4a9 
i : 2 
ao 5) ay. 2 4ao (4.103) 


From Section 4.3 using equation (4.101) 


dN(x) 2 
aga ce ex 
OL =a, 
B = 2a,f-a, f as large as needed (4.104) 


These results give two bounds, one for a positive ¢ and another for a 


negative €. 


f-a, y)+ 
2 yee I ll < | L 2 | es 
a 
ex<O — |x || < [i521] (4.105) 


The results obtained by Barrett were also obtained by Lepschy, 
Marchesini and Picci [17] using a comparison method while a slightly 
smaller region was obtained by Christensen [4], using contraction mapping 


and fixed point theorems. 
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CHAPTER V 


INITIAL CONDITIONS AND THE VOLTERRA SERIES 


Dot Introduction 
All of the examples with the Volterra series used up to here 
have assumed zero initial conditions. In this chapter the inclusion 
of initial conditions of the differential equation describing the 
linear plant will be discussed and a method for incorporating the initial 
conditions into the series will be developed. This method, coupled with 
the previous method for evaluating the series, gives a relatively simple 
method of determining the Volterra series for non zero initial conditions. 
Using the initial condition method, it is then shown how to 
treat a step input to give a quickly converging series. In practice, 
any input which will tend to a non zero average value can be treated 


in this way. 


5.2 Inclusion of Initial Conditions in the Volterra Series [5] 


Consider the general system in Fig. 5.1. 





Fig. 5.1 The General System 
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where L and L are linear operators, and N is a single-valued time- 
invariant nonlinearity. Let 


c sh l4. te 


L (s) = GS-1) 
sP+a Bat ea 
1 P 
dst "4+.. .4d 
LCs 3 oF (592) 
st4b rat ae 
1 q 
_ i 
N(x) = } n, x) (5.3) 
i=2 


and let L. and Ly have initial conditions TeCOLyaee WOL ye eee 2 (OL), 
and eG ee Cees ye Ory} respectively, where y(0+) is the 
value of y(t) just after t=0, and yD (t) is the q'th derivative of y(t). 


Consider first the equation 
yi a= L, (8) (524) 


The solution to this equation, including initial conditions, is given in 


Laplace transforms by, [8] 


4 es : ey S-)) (04) 
al cms GaAs] 
AR) oo ee ee ee L(s)G(s) (5.5) 


Beh an | gus 
fh. q 


where the following definitions are made. 


by = 1 


y © (o+) = y(0+) (5.6) 
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Equation (5.5) can be written in the time domain as 


y(t) = yo(t) + Ly (e) | (5.7) 


where y, 6) is the initial condition function given by the inverse 
transform of the appropriate function in equation (5.5). 


In a similar fashion the equation 


(tees L, Ce) (5.8) 


can be treated 


Dolrep& 


Tawed (acd) 
A s esas x (O+) 
+ L,(s)E(s) (5.9) 


X(s) 
where similarly 


a. = 1 and x 9) (ory = x(0+) | (5.10) 


Taking the inverse Laplace transform gives 


CC) oer CL) eer L, fe) (5.11) 


Now expressing the system output as x(t) and the input r(t) gives 


MG) oe ate OS Re eee = L,#N(x)) (5-12) 
saaxeit)s + L(r(t) - y,(t))- LN (x) (5.13) 

where as before 
L twieL We (5.14) 
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Equation (5.13) is the equation representing the general system 
under consideration with initial conditions on the linear plants. The 
effect of initial conditions is to modify the input function. Clearly all 
the previous results obtained about stability of the system and the 
convergence of the Volterra series will also hold in this case for the 
modified input function. 

In the evaluation of the Volterra series the modified input 


function could be expressed as 
= ‘ + = 
w(t) = x, (t) +L (r(t) - ¥,(0)) (5.15) 
and then the output could be expressed as a function of this new function 


we CG) « 


%(t) (5.16) 


i] 
< 
= 
ct 
ww 
! 
{7 
* 
{2 
= 
tal 
wv 


1] 
ie | 


x(t) H, (w) AG) 


i=l 


This appears to simplify the kernels of the series a little but 
i 
in effect just removes the factor TT tits) from the expression for each 
k=1 


H, (sy>+++S4)- 


The case where N(x) contains a linear term can be treated the 


same way as was suggested in Section 4.5. 


5.3 Example of Initial Conditions 


Consider the system of Fig. SL with 


tats) “355 (5.18) 
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Equation (5.27) now becomes 


Ce) ge LY ge KES Ge 


and expressing x(t) as a series function of w'(t) gives 


x(t) = H(w') = } H,(w') (Ss 
j=l ° 


As before, substituting equation (5.32) into equation (5,31) 
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5.4 Modifications of the Series for Step Inputs. 


Consider the example worked out in Section 3.6 of a system 
with a step input. In this example it was found that with inputs of 
a small magnitude there was good convergence of the series for all 
time. However, as the input became larger the series converged for 
a finite time interval but diverged quite rapidly after that. It is 
known from computer studies that this system is quite stable for all 
inputs and so the question is raised of whether it is possible to 


transform the equation so that the series converges for all inputs. 
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(5.40) 


A method for achieving this is outlined below which uses a transformation 


and the initial value method of the previous sections. 


The system to be considered is that of Fig. 5.1. 


(Ee mre) tat eg (48 een (t)) eee NG) (5.41) 


with L_, L, and N as defined in equations (3... 20 5n2) 5 (228) 
<a 1) re 


x (t) and y_(t) are initial condition functions as defined in Section 
oO oO 


5.2 and 


L = L +L (5.42) 
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Assume both linear operators are stable, so that the initial condition 


functions tend to zero as t increases. 


lim x(t) 


[Eeae8 


It 
io) 


i} 
So 


lim y(t) 


exaee 


lim L(y, (t)) = 0 (5.43) 


treo 


Also assume that for an input function 
r(t) = Uu(t) (5.44) 


that the output is bounded, and that the output reaches a steady state 


value. Then, 


IIx] < Ky 
a [Nx] < K, (5.45) 
lim x(t) = X- (5.46) 
Ge 
Let 
L(x) = (fo h(t )x(t-t) det (5.47) 


and from the assumptions of stability of the linear operators, then 
|n(t) | < K,exp(-vt) y>0 (5.48) 
Consider the limit as t>~ of L&N (x). 


LAN(x) = [© h(T)NGe(t-1)) dt (5.49) 


aoktythaos [aistok 4d 









S : 
ome ajasyeteqo room “died 9m ae ap 


Wh as 
a 


‘ : : aon 
-29ek9 me” 4 ah os 09 hiss enotson 7 
, 7 - 1 


0 = €3) ox ot! gees | 
o a ; 
} _ 7 
an 
0 =- (3) _y tl 
G 
wd} - 
Ua 
((9) ¥)_ J mil a 
o°"s of 
. oa 
juqni me toi jedgd omen oalA 
(ayn = - 
‘odd gadt base ,bebaved @& apgauo ens 9 
- 5 gi 





ee, 
nett nates a 





- 


gE 


5 ee 


70 
It is assumed that both h(t) and x(t) are equal to zero for all t<0, so 


that equation (5.49) may be written 


LAN(x) = [0 h(t)N(x(t-t))dt (5.50) 
and then 
lim L*N(x) = als h(t) (lim N(x(t-t))) dt (5.51) 
t?-o too 


since both |n¢t) | and |N (x) | are bounded. N(x) is single valued and 


continuous, so that using equation (5.46), equation (5.51) becomes 


lim L*N(x) 


treo 


of RCA)NG,) dt 


i} 


Nee h(t)dt Cy a52) 


Now AP h(t)dt exists, and so using Laplace transform theory 


[8] to evaluate the integral gives, 


ie h(t)dt lim 4 h(t)dt 


t+ 


= lim-sL(s)/s 
so 


=_lim L(s) (5.353) 
s>o 
Incorporating this result into equation (5.41), and taking-the Limit 


as t>~ on both sides, gives 
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If the algebraic equation (5.54) has a real solution for Xe then 
this is the final value of x(t). So now in equation (5.41) substitute for 
x a function z(t) + Xe Then z(t) is a function whose steady state value 


is zero and has an initial value given by 


z (0+) x(0+) - 5 nein ole (5255) 


z (0+) x(O+) i oS 


and so on for the initial values of the higher derivatives. Making 


these substitutions gives 


PACE Se GS z(t) + U ae Ls) - Ly, (t) - L*¥N(z(t) + x) 


(5.57) 


z(t) - L, y,6t) + U lim L,(s) - N(x) lim L(s) - 
s*o s>o 


L¥N' (2) 


=" z(t) a(t)) — eas CL) LAN (@) (5.58) 


where z(t) is the initial condition for x(t) with x(0+) replaced by 


z (O+) + x, and N'(z) can be written as 
M al i ; ee 
NU a toy Cue xe (5.59) 
= P dg j f£ 
i=2 j=l 


and where 


Cc. = 41! Giaap7 (5.60) 
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This new nonlinearity contains a linear term which can be 


removed and placed with the linear plant. 


4 al i=] 
z(t) = Zee) Layee) Le } ny Gp Bee 
i=2 
M aL i : 
eaUee) eer ce Meare te (Snel) 
f=2)5j=20 °° 
Tee 
I+ i = R ere 
(L on ix, Ly) 2X) z(t) Ly, 6h) L*N"(z) (5262) 
= -1 “1 el 7 
oe z(t) = K 2i(t) - K *L y(t) - K *LAN"(z) (5.63) 
where 
M a ij 4 aM 2M : ae 
N"(z)= } ni, eG. 2: Xp ae. ci elt eC! ey ne (S 664) 
i=2 j=2 J 12 jai ingens 
a : 1-1 
iS co aulisay Gamer as ome (5.65) 
i=2 


and x ley, is now the new linear system. It can be seen from the 
expression for N'(z) that the greatest power of the nonlinearity has not 
increased but more terms have appeared at lower powers. 

z(t) can be calculated from equation (5.63) in the form of a 


Volterra series. 


co 


z(t) = } oe (z(t) Leen) t) (5.66) 
n=1 


and then x(t) can be found from 


(2) 


sid honk gp ba CEE Poems we H (z(t) - Ly (t) ) (5.67) 
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and it can be seen that Xo plays the role of a zero order operator. 


5.5 Example 
Consider the same example as in Section 3.6. 


The system was 


NGO) = ex? 

1 
£6) ~ Sha 
L,(s) = 1 
x(O+) = 0 
w(t) = kit) 


Equation (5.41) becomes 


EE) @= Be) = ek) 
Taking the limit as t>~., 
3 
Xe = (k/a) -(¢/a) (x,) 
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From the solution of cubic equations, it is known that 
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which is the only real solution as long as 
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For stability of the linear system a>o and it then depends on 


conditions (5.73) to hold. In fact all that is needed is 


a> De 27K 
Now let 
xt) = 72 (Ct) + x, 


and substituting in equation (5.69) and cancelling terms 


ae ee ee 3) 2 2 
2(t)+x. = z(t)+C) (7) 5 eL (2 +3x,2 +3X¢ Zz) 
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z(t) z(t) - 3ExX, L, (2) - cL eC + 3x_2°) 


z(t) 


For this case using equation (5.9) 
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(5.74) 
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For convenience in notation, let 


(5.83) 


y = Semeera <1 (5.84) 


Then solving equation (5.77) for z(t) in the form of a Volterra series 


z(t) =) 2, (t) (5.85) 
jal 7 
2, (t)= w(t) = - eee 
25 (t)= -x(y/6) Gao eo + ae 5) 
24(t)= =x (92172) (472 **-1206 ***+102e° 326“ T43e ) 
2, (t)= mx (43/2160) (1069e" "43206 468552" 
536007 421156 3840 T4252 (9°) (5.86) 


These results were obtained by methods outlined before but grouping 
terms with the same power of x, or y under the same order of operator. The 
results are plotted in Fig. 5.2 where it can be seen that there is good 


convergence of the series. Compared to Fig. 3.1 there is more error for t 
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small but there is good convergence here at all inputs while the 
solution shown in Fig.3.1 does not converge for all time t and all 


inputs. 
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CHAPTER VI 


COMMENTS AND CONCLUSIONS 


6.1 The Volterra Series 

The Volterra series is a very useful tool for the nonlinear 
systems engineer. It can be used for a large class of systems whose 
nonlinearity can represented with sufficient accuracy as a finite 
power series, and some derivatives of this power series, when these 
derivatives can be incorporated as zeros of the linear plant. This is 
quite a large class of systems, as physical nonlinearities are usually 
not too abrupt, and input signals to systems are usually limited in 
magnitude by physical considerations. 

The main problem with the series, as with all types of power 
series solutions, is to obtain convergence of the series. For small 
inputs and the types of systems considered here, the series always 
converges but as the signals increase in magnitude, the series begins to 
diverge as is shown in the example in Chapter III. This does not mean 
that the Serene not stable however, and if it can be shown that the 
system isin fact stable then the system may be able to be transformed 
to a new equivalent system whose Volterra series does converge. This 
has been done for the case of the step input in Chapter V. 

Another method of evaluating the Volterra series by differen- 
tiating the series and substituting into the differential equation has 
been used by Flake [10], which should be able to handle all the types of 
nonlinearities that contain sums of powers and products of derivaties. 


A special method for including initial conditions is also given but there 
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is more work involved in calculating the series for these methods. 

Using the kernels of the series, it is possible to analyse the 
response of the system to quite a few inputs which engineers are interest- 
ed in. For sinusoidal inputs the distortion terms are quite readily cal- 
culated substituting s = tyw into the transforms of the kernels [11]. For 
initial conditions or step inputs then the methods developed here give 
good convergence of the series and for impulses the kernels themselves 
give the output function. 

One way of synthesizing a nonlinear system would be to use a 
general nonlinearity and work out several terms of the series. Then the 
coefficients of the powers of the nonlinearity could be chosen to give 


a desired result or close to it. 


6.2 Stability Analysis 


The methods of analysis used here to investigate stability 
depend on the functional analysis approach using contraction mapping. 
This approach combines well with the Volterra series and gives useful 
results. For instance, over the region in which the series converges 
absolutely, if the initial condition function is an exponentially 
damped fonctions chen by consideration of the series it can be seen 
that the output will also be exponentially damped and so the system 
is asymptotically stable [5]. 

For the expanded region of convergence as found here there is 
no guarantee of asymptotic stability in general. It appears from other 
work (see [1]) that all systems of first and second order with non- 
linearities as used here which are stable at the origin are asymptotically 


stable over the whole region of bounded-input bounded-output stability. 
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The problems arise with higher order systems with zeros in the linear 
plant. Lt Proutd be possible to determine the asymptotic stability 
bounds using the functional approach with or without the Volterra series 


and more work is needed in this area. 


6.3 Conclusions 

In this thesis a class of nonlinear systems is considered. 

It consists of up to two linear plants and one single-valued einte 
invariant nonlinearity which can be represented to a sufficient degree 
of accuracy by a finite power series. With this class of systems the 
easiest way to calculate the series was demonstrated with and without 
initial conditions and a transformation introduced to guarantee con- 
vergence of a stable system with a step input. 

The bounded-input bounded-output stability of this class of 
systems was investigated and for first and second order systems a 
sufficient condition was found that the slope of the nonlinearity eusaid 
be contained in the Hurwitz sector of the system with the nonlinearity 
replaced by a linear gain. For higher order systems, a bound was found 
on the output from which stability could not be assumed, but it has been 
found from computed results that for nonlinearities with slopes within 
the Hurwitz sector, the systems were bounded-input bounded-output 
stable so that more investigation is needed in this area. 

Also some work is needed on the prediction of limit cycles 
along with the bounded-input bounded-output stability. If it.can be 
shown that a system is bounded-input bounded-output stable and no limit 
cycle can exist then the system must also be asymptotically stable. This 


is especially important for the higher order systems. 
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In summary this thesis shows the use to which the Volterra 
series can be put in the analysis of nonlinear systems. The Volterra 
series was shown to be unique so that it can be calculated the way that 
is easiest and a proposed method of calculation was set out. A simple 
way of including initial conditions in the series is introduced and it 
is used to give faster convergence for step input functions. The con= 
traction region has been expanded for a large class of systems to give 
a convergent series for a finite time interval. Then using a diagonal- 
ization technique on first and second order systems bounded-input bounded- 
output stability is obtained for a nonlinearity whose slope lies in the 
Hurwitz sector of the system. Although this has been done by other 


methods, it is now possible to find a solution to the system. 
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